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Abstract 

Using ab initio calculation, we investigate systematically the structural and electronic prop- 
erties of N^Nbi+^Sni-.,; (x — 0, 0.25, 0.50). Here, projector augmented wave approach (PAW) 
implemented in the Vienna ab initio simulation package (VASP) within generalized gradient approx- 
imation (GGA) for the exchange-correlation functional has been used. In this article, it is reported 
that though Ni2NbSn and Ni2Nb1.25S1io.75 have no structural transformation, NiaNbi.sSno.s can 
transform to tetragonal structure from cubic L2i phase. The cubic lattice parameter decreases 
with Nb doping at Sn sites in off-stoichiometric alloys. The alloys are in paramagnetic phase in 
all the structures. The hybridization between Ni and Nb 3d states triggers the tetragonal distor- 
tion. Due to Nb doping in cubic L2i phase, there is a significant change in total density of states 
(DOS) at Fermi energy (Ep) (N(E^)). N(Ej?) increases with increasing Nb doping. But, N(E^) 
decreases during structural transformation of Ni2Nb1.5Sno.5- The superconducting critical temper- 
ature (Tc) also changes with Nb doping in cubic phase and tetragonal distortion because Tc very 
much depends on N(Ejr). 



1 INTRODUCTION 



Heusler alloys are ternary intermetallic compounds which can be defined by stoichiometric formula 
X2YZ, where X and Y are transition elements and Z is a III, IV, or V group element [I]. Cu2MnSn was 
the first Heusler alloy discovered by F. Heusler in 1903 [2]. Subsequently, a large number of Heusler 
alloys have been discovered with a wide range of physical properties. These class of materials present 
different magnetic properties like itinerant and localized magnetism, antiferromagnetism, helimag- 
netism, Pauli paramagnetism or heavy- fermionic behavior [HI [U El El El El El EU1 EH E2 E3 El] ■ Along 
with magnetic diversity, the alloys have acquired recent interest depending on electrical transport 
properties [6], shape memory effect [15) . half- metallic ferromagnetism [16], and semimetallic [H E] 
behavior. Moreover, several Heusler compounds with 27 valence electrons have been discovered to 
have a superconducting ground state [171 EHJ EU EOJ Ell E21 El]. Recently, T. Klimczuk et al. has 
listed wide varity of superconducting Heusler alloys [24] . 

One of the superconducting Heusler alloys is N^NbSn. This is crystalized in cubic L2i phase [23] in 
room temperature which is same as the other Ni based Heusler alloys [25] , Previously, some research 
groups [231 ESI ET] have mentioned the superconductivity, specific heat and electrical properties of 
Ni2NbSn. Among all the reported Ni based superconducting Heusler alloys [23], this stoichiometric 
alloy has highest Tc of 3.4 K [SUEZ]- Unfortunately, there is no reported data available about the 
properties of off-stoichiometric Ni2Nbi+ x Sni_ x alloys. So, the alloys have a lot of research scope in 
experimental as well as theoretical point of view. The aim of the present paper is to understand the 
structural and electronic properties of N^NbSn and its off-stoichimetric alloys theoretically. There are 
some fundamental questions, whether Ni-Nb-Sn alloys have structural and magnetic transitions and 
how the Tc varies in off-stoichiometric alloys. So, here the theoretical ab-initio calculation is adopted 
to find out the change in structural and magnetic properties as well as the variation of DOS of atomic 
orbitals with Nb doping at Sn sites of off-stoichiometric alloys, because change of value of the DOS at 
Ep changes the Tc appreciably. 

2 COMPUTIONAL DETAILS 

In the last two decades, tremendous progress in the development of methods for ab initio calculations 
have been witnessed to understand electronic structure and material properties. This development 
was based on density- functional theory (DFT) , which considers the complexity of the electron-electron 
interactions in many-electron systems into an effective one-electron potential, which is a functional of 
electron density only [28], [29] . 

One of the most accurate calculation can be done by VASP, developed by Georg Kresse and his 
coworkers [301131) . In this paper, the electronic structure of Ni-Nb-Sn alloys have been investigayed by 
the VASP. VASP is a plane-wave code for ab-initio density-functional calculations. It tries to match 
the accuracy of the most advanced all electron codes by using PAW approach (32^ [33] for describing the 
electron-ion interaction. Furthermore, the electronic exchange and correlation energies are described 
by using the GGA of Perdew et al. JHU [35]. In VASP, pseudopotentials have been introduced to avoid 
the need for an explicit treatment of the strongly bound and chemically inert core electrons. 

For electronic structure calculation of Ni2NbSn, cubic L2i structure having Fm3m space group 
with atomic positions 8c (Ni), 46 (Nb) and 4a (Sn) is used. Starting with the experimental lattice 
parameter a = 6.157^4 [23] of 16 atom unit cell, convergence test of KPOINTS and ENCUT has been 
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performed considering ground state energy difference (AEO) ~ 5 meV. The optimized plane wave 
cut off energy (ENCUT) of 450 eV and Monkhorst-Pack KPOINTS grid of 12x12x12 are employed 
for the calculation. Next, the lattice parameter (a) is relaxed so as to minimize the total ground 
state energy (E0)/unit cell. Afterward, Nb is doped at Sn sites in off-stoichiometric Ni2Nbi+ x Sni_a; 
according to x = 0.25 and 0.5. The lattice parameters of the off-stoichiometric alloys are also relaxed 
to minimize EO/unit cell. Using the equilibrium structure, spin integrated total DOS per unit cell 
and partial density of states (PDOS) per atom are calculated. The input cubic L2i crystal structures 
of Ni2Nbi +x Sni_ a; for x = 0, 0.25 and 0.5 are shown in Fig. [TJ Here Nbl and Nb2 indicate the Nb 
atoms which are situated at the parent Nb sites and at the Sn sites respectively. 

In Fig. [JJ the particular doping combination of Nb at Sn sites of Ni2Nb1.25Sno.75 and Ni2Nbi.sSno.5 
are equivalent to other possible combinations. Results remain unchanged for all possible configurations. 
So, the present doping combinations are reliable to use for calculations. 

3 RESULTS AND DISCUSSIONS 

3.1 Structural and magnetic properties 

First, we shall concentrate on the cubic structure of Ni2Nbi+ x Sni_ x . In Fig. [2] ground state energy 
E0 (eV)/unit cell of Ni2NbSn, Ni2Nb1.25Sno.75 and Ni2Nbi.5Sno.5 with respect to volume of cubic L2i 
structure are plotted. The theoretical EO/unit cell vs. volume curve is fitted with Murnaghan equation 
of state [36j [37], [38j [39] . It is observed that the equilibrium volume as well as the theoretical lattice 
parameter of Ni2Nbi +:r Sni_. r decrease with doped Nb concentration (x). 

The equilibrium unit cell volume and theoretical lattice parameters of Ni2Nbi +:r Sni_, r for x = 0, 
0.25 and 0.5 are listed in the Table [TJ 

From Table [TJ it is observed that the equilibrium lattice parameter of unit cell of Ni2NbSn is in 
good agreement with experimental value [23] within 1 %. Unfortunately, there is no experimental 
data available for lattice parameters of off-stoichiometric Ni-Nb-Sn alloys. It is obtained that, the 
theoretical lattice parameter of unit cell of Ni2Nbi +x Sni_ :r decreases with increasing concentration of 
Nb doping at Sn sites. Due to smaller size of Nb than Sn, the unit cell volume of off-stoichiometric 
alloys become less than the stoichiometric one. The change of theoretical lattice parameters with doped 
Nb concentration is depicted in Fig. [3j It is observed that the lattice parameter decreases linearly. 
So, from the trend of variation of theoretical lattice parameters, cubic unit cell lattice parameters of 
unknown compositions can be obtained. 

In the EO/unit cell vs. volume calculation, only the ionic positions of unit cell were relaxed for 
different cubic unit cell lattice parameters of Ni2Nbi+ a; Sni_ a: . Now, to find out whether Ni-Nb-Sn 
alloys have any structural transition, the ionic positions and volume of unit cells of Ni2Nbi+ a; Sni_. r 
( x = 0, 0.25, 0.5) are relaxed. Interestingly, it is obtained that Ni2NbSn and Ni2Nb1.25Sno.75 have 
no structural transformation. But, Ni2Nbi.5Sno.5 transforms to tetragonal structure where one axis 
elongates along one direction and the other two axes contract along the other two directions. For 
example, here a elongates and b & c contract. The equilibrium tetragonal structure of Ni2Nbi.sSno.5 
is of lattice parameters a = 7.31A, b = 5.70A and c = 5.70A Due to formation of Nb cluster in be 
plane the b and c axes contract and due to strong interaction between Ni & Nb a axis elongates. By 
this distortion, the system lowers its energy. In Table [2] the nearest neighbor distances between Ni & 
Nbl and Ni & Nb2 of Ni2Nbi.sSno.5 are listed in both cubic and tetragonal structures. 
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From the Table El it is observed that there are two different Ni-Nbl distances. The distance 2.64 
A of Ni-Nbl bond of cubic structure increases to 2.71 A due to tetragonal transformation. This Nbl 
is situated in the Nbl-Sn plane (see Fig. Q]). The distance 2.73 A of Ni-Nbl and Ni-Nb2 does not 
change due to structural transformation. These Nbl and Nb2 are placed in the plane of Nb cluster. 
The tetragonal distortion of Ni2Nbi.sSno.5 is quite similar to the ferromagnetic Ni2Mni.sSno.5 [40]. 
Here, Ni-Mn (at Sn sites) hybridization prompts the structural distortion. So, from the analogy, it 
is argued that hybridization between Ni and Nb (at Sn sites) states helps in tetragonal distortion of 
cubic L2i structure. 

In the following subsection, the total DOS and PDOS of Ni-Nb-Sn alloys in both cubic and tetrag- 
onal phases will be discussed to understand how atomic orbitals behave due to the Nb doping at Sn 
sites in cubic structure and structural transformation of Ni2Nb1.5Sno.5- 

Though Ni2Nbi +:r Sni_ a: for x = 0, 0.25 and 0.5 have such structural varieties, there are no magnetic 
transitions of these alloys. All the alloys are in paramagnetic ground state. So, magnetic moments 
are very less. This is contrast to the Ni based Heusler alloys. As elementary Ni is ferromagnetic in 
nature, Ni-containing Heusler compounds with high proportion of Ni are definitely expected to show 
magnetic order [4T]. But, here the presence of Nb states disrupts the magnetic order of Ni states 
giving rise to paramagnetism of these systems. Some research groups [23 |. 1271 14~2] have also mentioned 
the paramagnetic contribution to the magnetic susceptibility of Ni2NbSn. Waki et al. [23] identified 
experimentally that the magnetic susceptibility of these compounds shows temperature independent 
Pauli paramagnetism. According to them, the small amounts of magnetic susceptibility are consistent 
with the small values of DOS which is derived from specific heat measurements. From the values of 
DOS and temperature independent susceptibility, it is considered that the Fermi level is situated near 
a minimum in the DOS. Regrettably, there is no experimental data of paramagnetic susceptibility of 
off-stoichiometric Ni-Nb-Sn alloys in the literature. 

3.2 Density of states in cubic and tetragonal phase 

In Fig. |U El El H and El the spin integrated total DOS, PDOS of Ni, Nbl, Nb2 and Sn of Ni 2 Nbi +a; Sni_ :z; 
(x = 0, 0.25 and 0.5) in cubic L2i structure are depicted respectively. 

From Fig. EJ it is obtained that the valence band width of total DOS is 11 eV. From the PDOS 
plots of different atoms, it is obvious that the contributions to total DOS from 4 eV to -3 eV are 
mainly due to Ni, Nbl and Nb2 ti g and e g states; but Ni s and Nbl and Nb2 s and p states are much 
smaller. From Fig. El and El it is observed that around Ep, above and below the Fermi surface the 
contribution to the total DOS is mainly due to Nbl and Nb2 3d states. According to J. Winterlink et 
al. [43], a saddle point at a certain high-symmetry point in the energy dispersion curve at E^ leads 
to a high DOS at Ep. These are referred to as Van Hove singularities (VHS) [33] . Superconducting 
Heusler compounds with 27 electrons exhibit saddle point close to Ep in the corresponding energy 
dispersion curve according to the van Hove scenario [35l |H1 ES] • Here, the VHS are situated around 
1 eV just above Ep. Those singularities are mainly contributed by Nb 3d states. This same situation 
arises in other superconducting Heusler alloys [42 1 143 1 146] . Furthermore, near Ep around -0.4 eV 
there exists a minima in DOS. In 1985, Waki et al. |23] predicted the same situation experimentally, 
which is consistent with paramagnetism of the system. From Fig. [3] and El it is observed that arond 
-1 eV energy, the states are mainly Ni e g states. In the energy range -1.7 to -3.1 eV, the states are 
contributed by Ni e g and t2 g , Nbl t<i g and Nb2 e g & t% g . Due to tetrahedron field of Nb atoms, 
degeneracy of Ni 3d orbitals splits into up-lying twofold e g states and low-lying threefold ti g states. 
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From Fig. HI El [6] and it is obtained that the states around -4 eV are composed of Ni and Nb s 
states. From Fig. [5] and El it is observed that the states around -6 eV and -10 eV, are contributed by 
Sn p and s states respectively. Moreover, the gap in the energy range -6.9 to -8.5 eV is the typical 
behavior of Heusler alloys [16] . 

Now, in this context to see the effect of Nb doping at Sn sites in the cubic structure of off- 
stoichiometric Ni2Nbi+ x Sni_ x , we shall concentrate on total DOS at and around Ep for x = 0, 0.25 
and 0.5. The values of N(E^) with x concentration are listed in Table [3j From inset of Fig. [3] and 
Table [3l it is clearly seen that N(E^) increases with Nb doping. This occurs because at E^, the total 
DOS is mainly contributed by Nb 3d states and Nb PDOS increases with Nb doping. The value of 
DOS at Ep is significant because the value of Tc is dependent on N(E^). According to BCS theory 

; where Qp, is the Debye temperature and U is the electron-electron 
attractive interaction which is constant for the present alloys. Debye temperature can be written as, 

@£> = jr^ (j^r^J 5 [58j . In this defination, if N is the number of primitive cells in any specimen of 
volume V, the total no. of acoustic phonon modes is N. In the Ni-Nb-Sn 16 atom unit cell, there is 

only one primitive cell inside. So, @d = ^' wnere * s ^ ne volume of the unit cell and v is 

the velocity of sound in the alloy which is constant for Ni-Nb-Sn alloys. From Table [H it has been 
observed that Vq decreases in off-stoichiometric alloys due to doping of Nb. As a result, Qp increases. 
So, owing to increase of both Qp and N(E^), the Tc increases in off-stoichiometric Ni-Nb-Sn alloys. 

Moreover, the VHS shift towards the higher energy side from ~Ep and the maxima of DOS decreases 
with increasing x concentration. Due to Nb doping, number of valence electrons increases in the 
system. So, doping with charge carriers does not only change the density in the system but also 
smears the density of states diminishing the singularities [49J . 

In Fig. [9]and[10l spin integrated total DOS and Ni, Nbl and Nb2 e g and t 2g PDOS of Ni 2 Nbi. 5 SnQ.5 
in cubic L2i and tetragonal structures are depicted respectively. To understand the effect of tetragonal 
distortion on DOS, here we shall concentrate on DOS around Ep. 

From inset of Fig. [9j it is observed that the N(Ejr) decreases due to tetragonal distortion. This 
happens due to the Nb 3d DOS at Ep decreases in tetragonal structure (see Fig. [T0|) . The values of 
N(E_p) of Ni2Nbi.5Sno.5 in cubic L2i and tetragonal structures are listed in Table HI 

Due to tetragonal transformation, the volume of the cubic unit cell remains almost unchanged. So, 
the Debye temperature becomes constant in both cubic and tetragonal unit cell structures. Finally, 
due to decrease of N(Ejr) from cubic to tetragonal phase, Tc decreases. 

From Fig. [9l it is also observed that, the Ni e g peak at -0.94 eV shifts to -0.44 eV due to 
tetragonal distortion from cubic phase. This happens owing to hybridization between Ni and Nb e g 
states. Moreover, it is observed that the VHS disappear due to tetragonal distortion. This happens 
because some occupied Ni e g states shifts towards higher energy side after structural distortion and 
become unoccupied. Increasing DOS above Ep causes the disappearance of VHS and the system 
becomes stable. 



03, T c = 1.130 o exp 



4 CONCLUSIONS 

In summary, we would like to emphasize the important results of this article more precisely. Firstly, the 
theoretical cubic unit cell lattice parameter of Ni2NbSn is in well agreement with experimental value 
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within 1 %. So, the theoretical lattice parameter values of off-stoichiometric Ni2Nbi+ x Srii_ x is ex- 
pected to be also good agreement with the experiment. Second important point is that due to increas- 
ing concentration of Nb doping at Sn sites the theoretical cubic lattice parameters of off-stoichiometric 
alloys decrease. Due to smaller size of Nb than Sn, the unit cell volume of off-stoichiometric alloys 
become less than the stoichiometric N^NbSn. The trend of decrease of lattice parameters is linear 
in nature. So, from this tendency of linear behavior of lattice parameters, cubic unit cell lattice 
parameters of unknown compositions can be obtained. Next valuable result is that Ni2NbSn and 
Ni2Nb1.25Sno.75 have no structural transformation if the cubic structure is allowed to be fully relaxed. 
Interestingly, Ni2Nbi.5Sno.5 transforms to tetragonal structure where one axis elongates along one 
direction and the other two axes contract along the other two directions. Furthermore, it is observed 
that there is no magnetic phase transition along with structural transformation of Ni-Nb-Sn alloys. 
All the alloys are in paramagnetic ground state. The presence of Nb atom interrupts the ferromagnetic 
order of Ni atom. This gives rise to paramagnetism of the system. These results are very important 
because till now no experimental and theoretical information of structural and magnetic properties of 
off-stoichiometric Ni2Nbi +a; Sni_. r have been discussed in literature. Moreover, density of states plots 
of Ni2Nbi+ :r Sni__ a ; for x = 0, 0.25 and 0.50 give the idea how the electronic structure changes with Nb 
doping. First interesting result is obtained that VHS just above ~Ep, shift towards the higher energy 
side from ~Ep and the maximum value of DOS decreases with increasing Nb doping concentration. 
The VHS are mainly contributed by Nb 3d states. Doping with Nb increases the number of valence 
electrons in the system and increasing number of charge carriers smears the DOS out of Fermi sphere 
eliminating the singularities. It is also obtained that N(Ejr) increases with increasing concentration 
of doped Nb. This happens because at ~Ep, the total DOS is contributed by Nb 3d states and Nb 
PDOS increases with Nb doping. Now, Tc very much depends on N(Ejr). BCS theory [17] predicted, 



Tc = 1.130/j exp 



Due to increase of both ©£> and N(E^) with increasing concentration of 



UN(Ep) 

Nb doping, Tc increases in off-stoichiometric Ni-Nb-Sn alloys. This result is also completely new and 
increasing Tc has a great interest in making high temperature superconducting materials. Further- 
more, due to tetragonal distortion of Ni2Nb1.sSno.5j N(E^) decreases as Nb 3d PDOS decreases. In 
this context, one important point should be mentioned that hybridization between Ni and Nb 3d states 
triggers the tetragonal distortion of Ni2Nb1.5Sno.5- Due to tetragonal transformation, the volume of 
the cubic unit cell remains almost unaltered. So, the Debye temperature becomes constant in both 
cubic and tetragonal unit cell structures. Finally, due to decrease of N(Ejr) from cubic to tetragonal 
phase, Tc decreases. Moreover, the VHS disappear due to tetragonal distortion. This occurs because 
some occupied Ni e g states shifts towards higher energy side after tetragonal distortion and become 
unoccupied. Increasing DOS above F>p causes the diminish of VHS. Due to this mechanism, the dis- 
torted structure becomes stable of Ni2Nbi.sSno.s. Thus, this paper gives a deep insight of structural 
and electronic properties of superconducting Ni2Nbi+ x Sni_a; Heusler alloy. 
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6 Tables 



Table 1: List of equilibrium unit cell volume, theoretical lattice parameters of Ni2Nbi+ x Sni_ a; for x 
= 0, 0.25 and 0.5 in cubic L2i phase and experimental lattice parameter of Ni2NbSn. 



Doped Nb 
concentration 
(x) 


Equilibrium 
unit cell 
volume 

(i 3 ) 


Equilibrium 
theoretical 
lattice parameter 

(i) 


Experimental 
lattice 
parameter 

(i) 


0.00 


239.938 


6.214 


6.157 |23] 


0.25 


239.030 


6.206 




0.50 


238.201 


6.199 





Table 2: Nearest neighbor distances between Ni & Nbl and Ni & Nb2 of N^Nbi.sSno.s in both cubic 
L2i and tetragonal structures. 



Structure 


Ni-Nbl 

(A) 


Ni-Nb2 

(A) 


Cubic 


2.64, 2.73 


2.73 


Tetragonal 


2.71, 2.73 


2.73 



Table 3: List of theoretical values of N(Ejr) of Ni2Nbi_|_ x Sni_ x for x = 0, 0.25 and 0.5 in cubic L2i 
phase. 



Doped Nb 
concentration 

(x) 


N(E F ) 
(states/eV cell) 


0.00 


10.34 


0.25 


17.31 


0.50 


18.43 



Table 4: List of theoretical values of N(E^) of Ni2Nbi.5Sno.5 in cubic L2i and tetragonal phase. 



Crystal 
structure 


N(E_p) 
(states/eV cell) 


Cubic 


18.51 


Tetragonal 


14.99 
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7 Figures 




Figure 1: Cubic L2i structures of Ni2Nbi +I Sni__ :r for x = 0, 0.25 and 0.5. 
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Figure 2: Total ground state energy (EO)/ unit cell and fitted curve vs. volume plot of N^Nbi+^Sni-; 
for x = 0, 0.25 and 0.5 in cubic L2i structure. 



11 




Figure 3: Theoretical lattice parameter and fitted curve vs. doped Nb concentration (x) plot of 
r^Nbi+zSni-z for x = 0, 0.25 and 0.5 in cubic L2\ structure. 
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Figure 4: Spin integrated total DOS plot of N^Nbi+^Sni-^ for x = 0, 0.25 and 0.5 in cubic L2i 
structure. 
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Figure 5: Spin integrated Ni PDOS of N^Nbi+^Sni-,,; for x = 0, 0.25 and 0.5 in cubic L2i structure. 
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Figure 6: Spin integrated Nbl PDOS of Ni2Nbi +x Sni_a; for x = 0, 0.25 and 0.5 in cubic L2i structure. 
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;ure 7: Spin integrated Nb2 PD0S of Ni2Nbi +x Sni_ a; for x = 0, 0.25 and 0.5 in cubic L2i structure. 
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Figure 8: Spin integrated Sn PDOS of Ni2Nbi +:E Sni_ x for x = 0, 0.25 and 0.5 in cubic L2i structure. 



17 




-12 -10 -8 -6 -4 -2 2 4 



Energy (eV) 



Figure 9: Spin integrated total DOS plot of N^Nbi.sSno.s in cubic L2i and tetragonal structures. 
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Figure 10: Spin integrated Ni, Nbl and Nb2 3d PDOS plot of Ni 2 Nbi. 5 Sn .5 in cubic L2i and tetragonal 
structures. 
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